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Abstract 

A surface model with skeletons is investigated by using the canonical Monte Carlo 
simulations. The skeleton is composed of linear chains, which are joined to each 
other at the rigid junctions. A one-dimensional bending energy is defined on the 
linear chains, and no two-dimensional curvature energy is assumed on the surface. 
The model undergoes a first-order transition between the smooth phase and the 
crumpled phase. We conclude that the first-order transition of the surface model 
with skeletons is independent of whether the junctions are elastic or rigid. 
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1 Introduction 



A well-known surface model is the one of Helfrich, Polyakov and Kleinert 
(HPK) [1,2,3]. Tlie pliase structure of HPK model is connected to string mod- 
els and membranes [4,5,6,7,8,9,10] and therefore has long been investigated 
theoretically [11,12,13,14,15] and numerically [16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31]. 

Membranes are considered to be two-dimensional surfaces in the conventional 
HPK model, and hence homogeneity is assumed in the model. However, mem- 
branes are not always homogeneous; cell membranes are known to be sup- 
ported by skeletons. Hop diffusion of membrane proteins or lipids was exper- 
imentally observed recently [32]. The free diffusion of molecules is prohibited. 
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and the diffusion is localized in some domains on the cell membranes. The ori- 
gin of this is considered to be due to the cytoskeletons. Artificial membranes 
are considered to have skeletons because they are partly polymerized [33]. 

The phase structure of skeleton models should therefore be studied as a sta- 
tistical mechanical problem. Skeleton models for the cytoskeleton were al- 
ready investigated in [34]. A hard- wall and hard-core potential was assumed 
on the polymer chains with junctions, and the responses to some external 
stress and the compression modulus were obtained [34]. Giant fluid vesicles 
coated with skeletons was experimentally investigated, and the mechanical 
properties were reported [35], where the actin filaments introduce an inhomo- 
geneous structure in homogeneous artificial membranes. A compartmentalized 
surface model was reported to undergo a first-order transition [36], where the 
compartment boundary prevents the vertices from the free diffusion. In [37], 
the phase structure of a surface model with skeleton was investigated, and it 
was reported that the model has a first-order transition between the smooth 
phase and the crumpled phase. The interaction of the model in [37] is de- 
scribed by one-dimensional bending energy for linear chains (or bonds) and 
two-dimensional bending energy for junctions. The two-dimensional Gaussian 
bond potential is also assumed in the Hamiltonian. Thus, the interaction of 
the model in [37] seems a little bit more complicated than that of HPK model. 
One of the reasons of this is because the two-dimensional elasticity is assumed 
at the junctions. 

Therefore, it is interesting to see whether the transition of [37] occurs in more 
simplified skeleton models. One possible model is obtained from the model of 
[37] by replacing the elastic junction with a rigid junction. 

In this Letter, we study the rigid junction model by using the canonical Monte 
Carlo (MC) simulations and see how the transition of [37] occurs in such a 
simplified model. We must note that the rigid junction model is not included 
in the elastic junction models. In fact, there are two types of elasticity at the 
elastic junctions; one is the out-of-plane elasticity and the other is the in-plane 
elasticity. The former elasticity can be rigid in the elastic junction model of 
[37] in the limit of infinite bending rigidity bj ^oo, however, the latter one can 
not be controlled in the elastic junction model. Therefore, the rigid junction 
model in this Letter and the elastic junction model in [37] are considered to 
be two different models. 



2 Triangulated surfaces 

The model is defined on the triangulated surfaces, which are characterized by 
N the total number of vertices including the junctions, Ng the total number of 
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vertices on the chains, Nj the total number of junctions, and L the length of 
chains between junctions. The junctions are assumed as rigid plates; twelve of 
them are pentagon and the others are hexagon. It should be noted again that 
Nj is included in N; a junction is counted as a vertex. Figures 1(a) and 1(b) 
show surfaces of size {N, Ns, Nj, L) = (2322, 600, 42, 6) and (A^, A^^, Nj, L) = 
(6522, 1482,42, 11), respectively. Thick lines denote the chains terminated at 
the junctions, and the chains together with junctions form the compartment 
boundary. All the vertices can fluctuate only locally on the chains as well as 
inside the compartments, and they are prohibited from the diffusion because 
of the fixed connectivity nature of the lattice. 




(a) (2322, 600, 42, 6) (b) (6522, 1482, 42, 11) 

Fig. 1. Starting configurations of surfaces for MC of size (a) 
{N,Ns,Nj,L) = (2322,600,42,6) and (b) {N,Ns,Nj,L) = (6522,1482,42,11), 
where is the total number of vertices including the total number of junctions, Ns 
is the total number of vertices on the chains, Nj is the total number of junctions, 
and L is the length of chains between junctions. Thick lines denote the chains, 
which terminate at the junctions. 

The lattice of {N,Ns,Nj,L) = (2322,600,42,6) in Fig.l(a) corresponds to 
that of (A^, Ns, Nj, L) = (2562, 600, 42, 6) in Fig.l(a) of [37]. The reason why 
A^ = 2322 of the surface in Fig. 1(a) is smaller than A^ = 2562 of the one in [37] 
is because the junctions are rigid objects in this Letter while they are elastic 
ones in [37]. One hexagonal junction reduces A^ by 6, and one pentagonal 
junction also reduces A^ by 5 if they were assumed as rigid objects. Thus, we 
can check that 2562 = 2322 + 6 x 30 + 5 x 12, where 30 and 12 are the total 
number of pentagonal junction and that of hexagonal junctions, respectively. 

The triangulated surfaces such as the ones shown in Figs. 1(a) and 1(b) are 
constructed as follows: Firstly, we process the icosahedron by dividing every 
edge into ^-pieces of uniform length, and obtain a triangulated surface of size 
A^ = 10£^+2 (= the total number of vertices), where i is the number of division 
of an edge in the icosahedron. The obtained surfaces are thus characterized 
by = 12 and Nq = N — 12, where Ng is the total number of vertices with 
co-ordination number q. Those A^5(= 12)-vertices and A'6(= A^ — 12)-vertices 
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respectively form pentagonal junctions and hexagonal junctions together with 
their nearest neighbor vertices, in the next stage. Secondly, compartmentaHzed 
structures are obtained by dividing I further into m-pieces (m = 1,2, ■■•), 
and we have uniform chains of length L= (£/m) — 2. Those chains terminate 
at the above-described pentagonal or hexagonal junctions. These junctions 
include 6 (or 7) vertices and 10 (or 12) bonds according to whether they 
are pentagonal or hexagonal, and hence each junction has many degrees of 
freedom at this stage. Thirdly, the junctions are assumed as rigid objects on 
the triangulated surfaces: the pentagonal (hexagonal) junction is identified 
with a regular pentagonal (regular hexagonal) plate. As a consequence, the 6 
(or 7) vertices together with the 10 (or 12) bonds are identified with a pentagon 
(or a hexagon). Thus, we have compartmentalized surfaces with rigid junctions 
as those shown in Figs. 1(a) and 1(b). 

On the surfaces shown in Figs. 1(a) and 1(b), we have L = 6, 16 and L = 11, 
£ = 26, respectively. We also note that the compartmentalized structures shown 
in Figs. 1(a) and 1(b) are identical to those of the model in [36] except the 
fact that the junctions are rigid. The total number A^^; of the compartments 
depends on the surface size N . Consequently, Nc is increased with increasing 
N just as in the surfaces in [36] . We choose the chain length L between the 
junctions as constant so that it is independent of N in the model. We fix the 
chain length L such that 



which respectively correspond to the values n — 21, n = 66, the total number 
of vertices inside a compartment [36] . The reason why we fix n is because the 
size of compartment is considered to be finite in the cell membranes, and also 
because it is expected that total number of lipids in the compartment remains 
finite in the cell membranes. 



3 Model 

The Hamiltonian of the model is given by a linear combination of the two- 
dimensional Gaussian bond potential Si, the one-dimensional bending energy 
32-1 which are defined by 



L = 6, L = 



11, 



(1) 



Si = Y. {X, - X,f , S2 = E(l - cos %,■)). 



(2) 




In these expressions, in Si denotes the sum over all the bonds {ij) con- 
necting the vertices i and j, and in 'S'2 denotes the sum over bonds i and 
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J, which contain not only bonds in the chains but also virtual bonds that con- 
nect the center of a rigid junction and the neighboring vertices on the chains. 
The symbol 6'(ij) in 5*2 is the angle between the bonds i and j, which include 
the virtual bonds described just above. The Gaussian potential Si is defined 
not only on the chains but also on all other bonds. As a consequence, the 
model is considered to be a surface model, although the mechanical strength 
is maintained by one-dimensional elastic skeletons joined to each other at the 
rigid junctions. 

Figure 2 shows a hexagonal rigid junction connected to chains, where the 

angle ^(jj) is defined not only at the vertices on the chains but also at the 
corners {=virtual vertices) of the junction. The triangular lattices attached to 
the chains were eliminated from the schematic drawing in Fig. 2 for simplicity. 



Fig. 2. A hexagonal junction connected to chains. The angle in ^2 is defined 

not only at the vertices on the chains but also at the corners (=virtual vertices) of 
the junction. The triangular lattices attached to the chains were eliminated from 
the schematic drawing for simplicity. 

Wc must comment on the size of the junctions. The junctions are two-dimensional 
objects and therefore have their own size to be fixed. The size of junction can 
be specified by the edge length R; the perimeter length of the pentagonal 
(hexagonal) junction is therefore expressed by 5R {6R). In this Letter, we fix 
the size of the junctions such that 



The value R = 0.1 is quite smaller than that of the elastic junctions in [37], 
where the edge length squared is R^ ~ 0.5 because of the relation Si/N = 
1.5 satisfied in the model of [37]. Since the junctions are two-dimensional 
rigid objects, it is expected that large-sized junctions infiuence equilibrium 
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R = 0.1. 



(3) 
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properties of the surface. Therefore, in order to reduce such unclear size-effects 
we choose the size R relatively small such as in Eq.(3). 



We must note that the junctions in the snapshots of Figs. l(a),(b) were drawn 
larger than that expected from Eq.(3); the junction size in the snapshots is 
fixed according to 5'i/A'"~1.5. In the following section, we discuss how do we 
fix the size R to the assumed value in Eq.(3) in the MC simulations. 

The partition function Z of the model is defined by 



/ n dX, exp [-S{X)] , S{X) = 5i + 6^2, (4) 

i=l 



where b is the bending rigidity corresponding to the one-dimensional bending 
energy, and /' denotes that the center of the surface is fixed. The integration 
Hill dXi is a product of the integration over vertices and that of junctions 
such that 



N 

l[dX,= n dX, n dXi, (5) 

i=l vertices i junctions i 



where Iljunctions i dXi is the integration over the degrees of freedom for three- 
dimensional translations and rotations. 



The bending rigidity b has unit of kT, where k is the Boltzmann constant, 
and T is the temperature. The surface tension coefficient a of Si is fixed to 
a = 1; this is always possible because of the scale invariant property of the 
model. In fact, in the expression aSi + bS2 we immediately understand that 
a = 1 is possible, because the factor a of Si can be ehminated due to the scale 
invariance of the partition function. Since the unit of a is (1/length)^, the 
length unit of the model is given by y'T/a. We use the unit of length provided 
by ^Jl/a — 1 in this Letter, although a is arbitrarily chosen to be fixed. 

It must also be noted on the relation of the size R in Eq.(3) to the thermo- 
dynamic limit N^oo. As discussed in the last part of the previous section, 
the limit — oo is taken in our model so that the compartment size remains 
constant as in the model of [36] . As a consequence, R is also remained constant 
in the limit N^oo. 
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4 Monte Carlo technique 



The integration ni=i d^i in the partition function is done by the canonical 
Metropohs technique. The update of X in MC can be divided into two steps, 
which arc corresponding to the integrations Hj unctions i ^-'^i ^nd Hvertices i ^-'^i 
in Eq.(5). The first is the update of X for the vertices including those on 
the chains: X are shifted so that X' = X + 6X, where 6X is randomly cho- 
sen in a small sphere. The new position X' is accepted with the probability 
Min[l, exp(— A^")], where AS — 5'(new) — 5'(old). The second is the update 
of the position of the junctions as three-dimensional rigid objects. This can 
be further divided into two processes: the first is a random three-dimensional 
translation, and the second is a random three-dimensional rotation. All of 
these MC processes are independently performed under about 50% accep- 
tance rate, which is controlled by small numbers fixed at the beginning of the 
simulations. We introduce the lower bound 1 x 10~^ for the area of triangles. 
No lower bound is imposed on the bond length. 

The junction size R is fixed to i? = 0.1 in Eq.(3) during the thermalization 
MCS. The initial value of R is given by i?~0.7 on the surfaces such as those 
shown in Figs. 1(a) and 1(b). Thus, we reduce R from R — 0.7 to R = 0.1 by 
6 X 10"^ at every 25 MCS in the first 2.5 x 10^ MCS. Because of this forced 
reduction of the junction size, the equilibrium statistical mechanical condition 
seems to be violated in the first 2.5 x 10^ MCS. Therefore, relatively many 
thermalization (1.75 x 10^ or more) MCS is performed after the reduction. 
Further thermalization MCS should be performed, if necessary. 

We use surfaces of size (A^, Ng, Nj, L) = (5222, 1350, 92, 6), (9282, 2400, 162, 6), 
(14502, 3750, 252, 6), and (20882, 5400, 362, 6) for the length L = 6, and {N, Ns, Nj,L)^ 
(6522, 1200, 42, 11), (14672, 2700, 92, 11), and (26082, 4800, 162, 11) for the length 
L = 11. A random number sequence called Mersenne Twister [38] is used in 
the simulations. 



5 Results 



Snapshots of surfaces are shown in Figs.3(a)-3(d). Figure 3(a) is a surface of 
size (A^, Ns, Nj, L) = (26082, 4800, 162, 11) obtained in the crumpled phase at 
b = 12.3, and Fig. 3(b) is the one obtained in the smooth phase at b = 12.4. 
The surface sections of Figs. 3(a) and 3(b) are shown in Figs. 3(c) and 3(d), 
respectively. These figures are drawn in the same scale. We immediately see the 
surface in the smooth phase is actually swollen, while the surface is collapsed 
in the crumpled phase. 
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(a) Collapsed surface at 6=12.3 (b) Smooth surface at 6=12.4 




The surface section 



(d) The surface section 



Fig. 3. Snapshot of the surface of size {N,Ns,Nj,L) = (26082,4800,162,11) ob- 
tained in the crumpled phase at (a) 6=12.3 and in the smooth phase at (b) 6 = 12.4, 
both of which are close to the transition point. The figures are drawn in the same 
scale. 



Si/N 



1.518 



1.514 



(a) 




Sj/N 



1.512 



1.508 



(b) 



L=ll 



O :N=26082 
A:N=14672 
□ :N=6522 



12.5 , 13 

b 



Fig. 4. The Gaussian bond potential Si/N vs. b obtained on the surfaces of (a) L = 6 
and (b) L = ll. Si/N shghtly deviates from S'i/A^~1.5. The curves are drawn by 
the multihistogram reweighting technique. 

The Gaussian bond potential Si/N is shown against b in Figs. 4(a) and 4(b), 
which correspond to the conditions L = 6 and L = 11, respectively. The solid 
curves are drawn by the multihistogram reweighting technique [39] . The values 
of Si/N in the figures slightly deviate from S'i/iV= 1.5, which is satisfied on 



8 



the surface without the rigid junctions or with rigid junctions of neghgible size. 
The reason of this discrepancy is because the surface includes the rigid junc- 
tions of finite size. A vertex is the zero-dimensional point, while the rigid junc- 
tion is the two-dimensional plate and hence shares some area on the surface. 
We find a gap or a jump in Si/N of the {N, Ns, Nj, L) = (26082, 4800, 162, 11) 
surface in Fig. 4(b), which can be viewed as a sign of a discontinuous transition. 




100 

50 



L=ll 



N=26082 

N=14672 
i 



J 



/ N=6522 
iJs n 



(b) 



12.5 , 13 
b 



Fig. 5. The mean square size against b obtained on the surfaces of (a) L = 6 and 
(b) L=ll. The curves arc drawn by the multihistogram reweighting technique. 



The mean square size is defined by 



(6) 



where X is the center of the surface. The crumpling transition is conven- 
tionally understood as the one of surface fluctuations accompanied by surface 
collapsing phenomena, which can be seen in our surface model; we have seen a 
collapsed surface and a swollen surface in Figs.3(a)-3(d). Therefore, we expect 
that reflects the crumpling transition on spherical surfaces. Figures 5(a) 
and 5(b) are plots of X'^ against h obtained under L = 6 and L—11. We find 
that the variation of X^ becomes sharp against b as N increases. Thus, it is 
expected that the variation of X'^ has a jump at intermediate value of b in 
either case of L. 

The crumpling transition is originally understood as the one for surface fiuctu- 
ation phenomena. Therefore, the bending energy S2/Ng is expected to refiect 
the transition, where Ng is given by 

N's^Ns + 6Nj- 12. (7) 



S2/N'g is the bending energy per vertex, because N'g is the total number of 
vertices where S2 is defined. N'g includes virtual vertices which are the corners 
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of the junctions (see also Fig.2). Those virtual vertices are not counted as the 
vertices and hence are not included in A^^. Total number of virtual vertices 
are 6Nj — 12, because the hexagonal junction includes 6-virtual vertices, and 
the total number of pentagonal junction is 12. Thus, we have Eq.(7) for N'g, 
and therefore we have A^^ = 1890, A^^ = 3360, A^^ = 5250, and A^^ = 7560 
for the surfaces of length L = 6 and size (A^, Ng, Nj, L) = (5222, 1350, 92, 6), 
(9282,2400,162,6), (14502,3750,252,6), and (20882,5400,362,6); and A^^ = 
1440, A^^ = 3240, and A^^ = 5760 for the surfaces of length L = 11 and size 
(AT, Ns, Nj, L) = (6522, 1200, 42, 11), (14672, 2700, 92, 11), and (26082, 4800, 162, 11). 

Figures 6(a) and 6(b) are plots of S2/Ni, against b obtained under L = 6 and 
L= 11. We find the expected behavior in S2/Ng under both conditions L = 6 
and L = 11; S2/Ng has a gap (or a jump) at intermediate b. This clearly shows 
that the crumpling transition is of first order. 



S2/Ns' S2/Ns' 

— I — ' — I — I — I — ' — I — ' — I — I 




Fig. 6. The one-dimensional bending energy S2/Ng vs. b obtained on the surfaces 
of (a) L = 6 and (b) L = ll. 



The transition can also be reflected in the two-dimensional extrinsic curvature, 
which is defined by 

5f = E(l-n.-n,), (8) 



where rij is the unit normal vector of the triangle i. In 5*2 , Y,(ij) denotes the 
summation over all nearest neighbor triangles i and j that have the common 
bond (ij) , which includes bonds belonging to the skeleton chains. We denote 
the total number of bonds, where 82^^ is defined, by Nb, which is given by Nb — 

J2(ij) 1- Note that Nb includes virtual edges, which are the one-dimensional 

(2) 

edges of the rigid junctions. In fact, we define 5*2 on the virtual edges, because 

it is reasonable to define extrinsic curvature on those edges. It is also noted that 

(2) (2) 
S2 is not included in the Hamiltonian, and therefore S2 gives no mechanical 

strength to the surface. 
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Fig. 7. The two-dimensional bending energy ^2 /Nb against b obtained on the 

surfaces of (a) L = 6 and (b) L = ll. S2 is defined by Eq.(8) and is not included in 

(2) 

the Hamiltonian. Nb is the total number of bonds where S2 is defined. 

The two-dimensional extrinsic curvature S2 /Nb is plotted in Figs. 7(a) and 
7(b) against b, which are corresponding to the conditions L = 6 and L = ll. Nb 
is the total number of bonds described above. We find that the dependence 
of 5*2 /Nb on b shown in Fig. 7 is almost identical to that of in Fig. 6. 

The gap (or jump) seen in Sf' /Nb also supports that the transition is of first 
order. 



The specific heat corresponding to the one- dimensional bending energy S2 is 
defined by 



b 



2 



Cs, = ^{{S2-{S2)f), (9) 



which can also reflects phase transitions if it has an anomalous behavior. 
Figures 8(a) and 8(b) show against b obtained under L = 6 and L ~ 
11. Solid curves drawn in the figures were obtained by the multihistogram 
reweighting technique, and those curves clearly show an expected anomalous 
behavior indicating that C52 is divergent when Ng — > 00 (or equivalently 
00). 

The specific heat corresponding to the extrinsic curvature s!^^ in Eq.(8) can 
also be defined by 

which refiects the transition as C5, does. Curvature coefficient for C„(2) was 

(2) 

assumed to be 1, because ^2 is not included in the Hamiltonian and therefore 
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Fig. 8. The specific heat for ^2 against h obtained on the surfaces of (a) L = 6 
and (b) L = 11. is defined by Eq.(9). The error bars on the symbols are the 
statistical error, which is obtained by the binning analysis. 



Cs2(2) Cs2(2) 




Fig. 9. The specific heat C„(2) for S2 against h obtained on the surfaces of (a) L = 6 
and (b) L = ll. C^{2) is defined by Eq.(lO). The error bars on the symbols are the 
statistical error, which is obtained also by the binning analysis. 

the curvature coefficient is not. Figures 9(a) and 9(b) show C„(2) against h 
obtained under L = 6 and L — 11. We can see in C„(2) the same anomalous 

■^2 

behavior as in Cs2- 

In order to sec the anomalous behaviors in and C„(2) in more detail, we 

'-'2 

plot the peak values of them in Figs. 10(a) and 10(b) in log-log scales against 
Ng and N, respectively. The straight lines were drawn by fitting the data to 

C^r^iN'sp, C-(tr oc (7V)-^ (11) 

"^2 

where cti, (T2 are critical exponents. Largest three data were used in the fitting 
in the case L — 6. Thus, we have 
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Fig. 10. Log-log plots of (a) against N'g and (b) C™2f against N' obtained on 

the surfaces of L = 6 and L = 8. The straight lines arc drawn by fitting the largest 
three data of Cg"^ and C^(^f to Eq.(ll). The peak values and the statistical errors 

for the fittings were obtained by multihistogram reweighting. 

(7i =0.911 ±0.118, (72 = 0.608 ± 0.076, (L = 6), 
(7i = 1.51 ± 0.08, (72 = 0.948 ± 0.262, (L = 11). (12) 

(72 = 0.608(76) for L = 6 is inconsistent with the fact that the transition is of 
first-order, however, (7i = 0.911(118) is consistent with that. (7i = 1.51(8) and 
(72 = 0.948(262) under L—11 support the discontinuous transition. 



6 Summciry and Conclusion 

A surface model with skeletons has been investigated by the canonical Monte 
Carlo simulations. The skeletons are composed of one-dimensional linear chains 
and rigid junctions, whose size is chosen sufficiently small compared to the 
mean chain-length. The surface is a triangulated sphere and divided into a lot 
of compartmentalized domains, whose boundary corresponds to the skeletons. 
The mechanical strength of the surface is given by the skeletons. There is no 
two-dimensional curvature energy in the Hamiltonian, while one-dimensional 
bending energy is defined on the chains connected to each other at the junc- 
tions. The two-dimensional Gaussian bond potential is included in the Hamil- 
tonian just as in the standard surface model of Helfrich, Polyakov and Klein- 
ert. The surface model in this Letter is considered to be different from the one 
with elastic junctions, because the rigid junctions cannot be identified with 
the elastic junctions due to the property on the in-plane elasticity. 

The surface is characterized by (iV, Ns, Nj, L), which are respectively the total 
number of vertices including the junctions, the total number of vertices on 
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the chains, the total number of junctions, and the length of chains between 
junctions. The length of chains was fixed to L = 6 and L=ll, which correspond 
to n = 21 and n = 66 the total number of vertices in a compartment. 

We found that the surface undergoes a first-order crumpling transition between 
the smooth phase and the crumpled phase. The one-dimensional bending en- 
ergy S2 has a gap (or a jump) at intermediate bending rigidity b, and the 

(2) 

two-dimensional extrinsic curvature 5*2 also has a gap at that point. These 
imply that the surface fluctuations are considered to be a first-order transi- 
tion. Moreover, it is found that the mean square size also has a gap at the 
transition point. This implies that the surface-collapsing phenomenon can be 
viewed as a first-order transition. 

The results in this Letter together with those in [37] show that the first-order 
crumpling transition can be seen in the spherical surface model even when 
the mechanical strength is maintained only by skeletons, which are composed 
of hnear chains joined to each other at the junctions. Moreover, the order of 
transition is independent of whether the junction is elastic or rigid. 

It has been reported in [36] that the transition can be seen in a compartmen- 
talized surface model, whose mechanical strength is maintained by the two- 
dimensional curvature of the surface in contrast to the model in this Letter 
and that in [37]. Therefore, we can also conclude that the first-order transi- 
tion occurs in the compartmentalized surface model independently whether 
the surface is mechanically supported by the skeleton (= the compartment 
boundary) or by the surface. 

Fluidity of lateral diffusion of vertices can be considered by using dynami- 
cally triangulated MC technique as in [36]. Vertices freely diffuse inside each 
compartment on the fluid surfaces supported by the skeletons. It is interesting 
to see how fluidity influences the transition of the skeleton-supported model. 
Many interesting problems remain to be studied on the surface model with 
skeletons. 
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